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ABSTRACT 



j> ■ We present the first spatially resolved observations of a sample of 23 Mira 

stars simultaneously measured in the near-infrared J, H and K' bands. The tech- 
nique used was optical long baseline interferometry, and we present for each star 
visibility amplitude measurements as a function of wavelength. We also present 
characteristic sizes at each spectral band, obtained by fitting the measured visi- 
bilities to a simple uniform disk model. This approach reveals the general relation 
J diameter < H diameter < K' diameter. 

6 

Subject headings: stars: atmospheres, fundamental parameters, variables — in- 
strumentation: interferometers — techniques: high angular resolution 
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1. Introduction 

Mira stars are variable, pulsating stars in the asymptotic giant branch (AGB). They 
have pulsation periods of about 100 — 1000 days, during which their light curves exhibit 
large amplitude variations (as high as 10 magnitudes at visual wavelengths, becoming a 
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more modest ~ 1 mag in the near-infrared), and their spectral types vary by a few sub- 
classes. During the pulsation cycle, their photospheric radii are thought to vary by 5 — 50 %. 
Mira star atmospheres are geometrically extended and have very low effective temperatures 
(< 3000 K), hence they are a rich molecular composition. They exhibit high mass loss rates 
(a few 10 -6 M Q /yr), and the ejected material forms layers which become sites of circumstellar 
maser activity and dust formation. Thus, although the Mira evolutionary stage is relatively 
short (~ 5 x 10 5 yrs), it is a crucial one in order to understand a variety of phenomena such 
as high mass loss, the formation of planetary nebulae, and the chemical enrichment of the 
galaxy. 

Despite the existing large body of observational and theoretical work, many aspects 
of the physics of Mira stars remain poorly understood; in particular their detailed atmo- 
spheric stratification and composition, and their dominant mode of pulsation. This is to a 
large extent a consequence of their intrinsic complexity (treatment of 3D convection, pho- 
tospheric shocks, and rich photospheric chemistry), and most needed are new observational 
data that can help discriminate between the large number of theoretical models that have 
been proposed. 

Spatially resolved observations of Mira stars are clearly desirable, and can be easily 
obtained using high resolution techniques such as long baseline interferometry. The geo- 
metrically extended and chemically rich nature of Mira star atmospheres manifests itself 
in apparent physical sizes which depend on the wavelength of observation, especially when 
crossing regions in the spectrum where the opacity changes rapidly. This fact was first ob- 
served in the visible region of the spectrum, where sizes as much as 50% smaller are found in 
narrowband near-continuum filters than in filters contaminated by TiO absorption (Bonneau 
& Labeyrie 1973; Labeyrie et al. 1977; Bonneau et al. 1982; Haniff, Scholz & Tuthill 1995; 
Hofmann et al. 2000, 2001). 

Although originally expected to be better probes of the underlying stellar continuum, 
measurements at infrared wavelengths also exhibit non-trivial wavelength dependence. Syn- 
thesizing measurements obtained at different facilities, it appears that as a general trend, 
Mira sizes appear large in the visible (Tuthill et al. 1994; Haniff, Scholz & Tuthill 1995; 
Ireland et al. 2004), smaller in the near-infrared (Perrin et al. 1999; Hofmann et al. 2002; 
van Belle, Thompson & Creech-Eakman 2002, and references therein), and larger again at 
thermal infrared wavelengths (Mennesson et al. 2002; Weiner, Hale & Townes 2003a,b). The 
interpretation of this broad trend is that the outer layers (i.e. those layers above a level that 
could be called "the photosphere") are relatively obscure in the visible, relatively transparent 
in the near-infrared, and again relatively obscuring in the thermal infrared. The physical 
basis for this variation of opacity with wavelength can be modelled once the underlying em- 
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prical variation of diameters with wavelength is established by observations. The purpose of 
this paper is to show that there is such a variation through the near- infrared, as defined by 
the J, H and K' bands. 

Due to instrumental limitations, simultaneous measurements at different wavelengths 
are scarce in the literature, but are clearly the optimum way to explore the wavelength- 
dependent structure, since Mira stars are suspected of exhibiting important cycle-to-cycle 
non-repeatibility. Dramatic size increases (25 — 100%) have been observed between near 
(~ 1.1 — 2.3fim) and thermal (~ 3.1 — 3.8/zm) infrared measurements, both in narrow 
(Tuthill, Monnier & Danchi 2000; Tuthill et al. 2000) and broad bands (Mennesson et al. 
2002). Smaller but significant (1—20%) effects have been observed in narrow bands within the 
near-infrared K-band (~ 2.0 — 2.4/im), as the narrow filters preferentially sample portions of 
the spectrum corresponding to strong molecular bands or the adjacent continuum (Thompson 
2002; Thompson, Creech-Eakman & van Belle 2002; Perrin et al. 2004). Using the same 
portion of the spectrum as in this work, Weigelt et al. (2003) studied the visibility of the 
Mira star T Cep dispersed across the J/H/K bands, finding similar J and H sizes which are 
~ 10% smaller than the K-band size; as well as size differences within the K-band as large 
as 26%. 

Clearly, the richness of the size-wavelength relation holds the potential to be the best 
probes of atmospheric structure in these stars. In turn, as has been illustrated in the recent 
multi- wavelength modelling effort by Perrin et al. (2004), a complete model that includes the 
stellar photosphere, dominant opacity sources and possible additional emitting layers must 
be globally fit in order to also extract accurate stellar diameters, and derive conclusions 
about the pulsation mode. 

In this paper we present the first systematic survey of Mira star diameters measured 
near-simultaneously in the near-infrared J, H and K' bands. This new multi-color data set 
allows us in particular to derive apparent size ratios, which are immune to many sources of 
calibration error. In § 2 we describe the measurement technique and Mira sample observed; 
in § 3 we describe the data analysis procedure; in § 4 we present our results, both in the 
form of uniform diameter fits at each spectral band in order to illustrate the main trends in 
these wavelength regions, and as tables of calibrated visibility data for straightforward use 
in global fitting including other datasets; in § 5 we discuss our main results and places them 
in context of recent theoretical work. 
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2. Observations 

Observations were carried out at the Infrared-Optical Telescope Array (IOTA), a long 
baseline interferometer located on Mount Hopkins, AZ (see e.g. Traub 1998). The data 
presented here were obtained using the two telescopes available at the time, on two con- 
figurations which provide physical baseline lengths of approximately 21 m and 38 m, at 
geographic azimuth angles (measured East of North) of —4 deg and +18 deg respectively. 

The beam combination and fringe detection instrumentation is exactly as described in 
Millan-Gabet (1999) and Millan-Gabet et al. (2001). In short, afocal telescope beams were 
combined at a bulk optics beam combiner, and interferograms were produced by modulating 
the optical path difference (OPD) between the two beams by ±60/im. The fringe modula- 
tion frequencies are in the range 100 — 2000 Hz, chosen as an optimum trade-off between 
maximizing the photon signal and minimizing the readout time to freeze out atmospherically 
induced OPD variations during the fringe acquisition. The two outputs of the beam splitter 
were focused on two well-separated pixels of a NICMOS3 array, the readout of which is 
synchronized with the optical path modulation (Millan-Gabet et al. 1999). 

The NICMOS3 camera contains a set of three standard near-infrared J (Ao = 1.25/im, AA = 
0.28/zm), H (A = 1.65fim, AA = 0.30yum) and K 7 (Ao = 2.16fim, AA = 0.32/im) filters, where 
(Ao, AA) are the center wavelength and full-width at half-maximum respectively. Rapid 
switching between the three filters guaranteed near-simultaneous observations (within a few 
minutes) in the three bands. 

A typical observation consisted of 500 interferograms, obtained in a few minutes, fol- 
lowed by a measurement of the sky background. Target observations are interleaved with an 
identical sequence obtained on calibrator stars which are unresolved by the interferometer 
or have small diameters which can be estimated. This sequence was repeated for each of the 
three filters. 

The basic properties of the Mira stars observed are listed in Table 1, and a log of 
observations is given in Table 2. 

3. Data Analysis and Validation 

The data reduction and calibration procedures are also essentially as described in Millan- 
Gabet (1999) and Millan-Gabet et al. (2001). For each interferogram, a peak fringe visibility 
(the fringe amplitude relative to the mean flux detected) is estimated by (a) correcting the 
intrinsic detector non-linearity, (b) sky subtraction, (c) subtracting the normalized outputs 



- 5 - 



of the beam splitter, to enhance the SNR, and (d) fitting the central three fringes of the 
resulting interferograms in the time domain to a point source response template computed 
as the Fourier transform of the spectral bandpass used in the observation. 

As a result of the fringe fitting procedure, we obtain the following quantities for each 
interferogram: mean, fringe amplitude, position of the central fringe, and sampling rate. 
The data are flagged as being a false fringe identification if an individual fringe visibility is 
outside the range 0.0 — 1.0, or if either the central fringe position or sampling rate has a value 
which is more than three times the standard deviation obtained from the ensemble of the 
500 reduced interferograms in the observation. Finally, an entire observation is rejected if 
the standard deviation of the central fringe position values exceeds L/2, where L = 120 /zm 
is the OPD sweep for each interferogram. 

For each observation, the above procedure results in several hundred estimates of the 
fringe visibility, one per interferogram that survives the flagging procedure. A final estimate 
for the observation is formed as the mean of those values, with an internal statistical error 
given by the error in the mean. 

Calibrator stars are measured close in space and time to the target observations in 
order to estimate the instrument's system visibility. Calibrated visibility amplitudes are 
obtained by dividing the target visibilities by estimates of the system visibility at the time 
of the target observations formed by spline interpolation of the calibrator observations. A 
correction to the system visibility resulting from the finite but known size of the calibrator 
stars is applied. Standard error propagation is used to include in the calibrated visibility the 
error resulting from the division by the system visibility and the systematic error resulting 
from the uncertainty in the calibrator star sizes. 

The calibrator sizes and uncertainties are in Table 3; these values were estimated as 
follows: (1) if available, adopt a near-infrared direct measurement from the CHARM cat- 
alog (Richichi & Percheron 2002); or else (2) if available, adopt the value based on stellar 
atmosphere modelling of Borde et al. (2002) or Cohen et al. (1999); or else (3) if available, 
and if the star is a luminosity class I, II, or III, use the empirical size-color relations of 
van Belle (1999); or else (4) fit the available visual and infrared photometry with a stellar 
blackbody model, using the Fbol module of the getCal software 1 . 

In order to obtain characteristic sizes at each wavelength, the calibrated visibility data 
obtained at each spectral band were fit using a simple uniform disk (UD) model for the stellar 



1 getCal is a a software package for planning of interferometer observations produced by the Michelson 
Science Center at the California Institute of Technology (http: //msc . caltech. edu/sof tware/index.html) 
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brightness. Although it is well known that for Mira stars the UD is in most cases not a good 
representation of the true center-to-limb variation (CLV), reconstruction of the true CLV 
is not possible when only sparsely sampled visibility data are available. In near-continuum 
bandpasses however, as near-infrared filters are to a large extent, UD profiles are expected 
to approximate the CLV better than at other wavelengths (Scholz 2003). Moreover, all the 
visibility data presented here were obtained at spatial frequencies well inside the first lobe of 
the visibility function (except for one star, R Leo, as discussed in § 4), where departures from 
the UD intensity profile are virtually undistinguishable, and therefore baseline dependent 
errors are expected to be minimal in this dataset. 

In Figure 1, we show 4 representative examples (out of a total of 34 similar observations) 
of the data and fitting results. The examples are chosen to illustrate distinct aspects of the 
data. The top two panels illustrate the general trend found that J diameter < H diameter < 
K 7 diameter. The bottom left panel illustrates our ability to sample relatively low visibilities 
for bright sources. The bottom right panel shows an example where Earth rotation allowed 
sampling of a relatively large range of projected baseline length, illustrating the visibility 
data precision for observations within a single night; here the visibility fit RMS is 0.01, or 
3.8% of the (changing) model visibility, which is a typical value for the dataset (full range is 
1-10%). 

The night-to-night data consistency was evaluated by comparing the results obtained 
when the same star was observed on different nights, at similar pulsational phase and using 
the same baseline, in order to isolate the test from possible size variations and non-circular 
star shapes. There are 3 such cases in our dataset (x 3 bands for a total of 9 independent 
comparisons) and the result is that the night-to-night consistency in visibility is in all cases 
better than 8%, with mean values being 2.3%, 3.0% and 3.8% for J, H and K' respectively. 

Our data can also be compared with previous measurements from the literature. This 
comparison can only be made at K-band. Considering the most recent measurements (Perrin 
et al. 1999; van Belle, Thompson & Creech-Eakman 2002; Monnier et al. 2004) we find that 
(a) our diameters show no systematic positive or negative bias with respect to the literature 
values, and (b) the diameter values agree to 20% or better. Considering that no effort 
was made to compare data taken only at similar pulsational phases, this level of agreement 
appears reasonable, given for example that K-band Mira sizes can change by 10% for a phase 
difference of only 0.04 (Perrin et al. 1999). 
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4. Results 

Recognizing that potentially the best use of our J/H/K' measurements will be as part 
of a larger modelling effort which incorporates for each star all the available data across 
the spectrum, a goal beyond the scope of this paper, we provide in Table 4 the totality of 
our visibility measurements. Columns 1 — 4 give the star name, calendar /Julian observation 
epoch, and filters used. Columns 5 — 7 give the 2-dimensional spatial frequencies (u,v), and 
measured visibilities. 

In addition, in order to capture the essential features of the "size" —wavelength relation 
across the near-infrared, we provide in Table 5 the results of fitting our visibility data to UD 
models. Columns 1 — 4 give the star name, observation Julian date, the pulsation phase 2 
and an estimate of the spectral type for that pulsation phase (see § 5). Columns 5 — 7 give 
the J, H and K'-band best fit UD diameters. The fact that we have measurements at 3 
different wavelengths, allows us to make a relative measurement in terms of diameter ratios, 
given in columns 8 — 9, which are expected to be immune to some sources of systematic error 
(e.g. calibrator diameter uncertainty, seeing effects) and also represent the linear size ratios 
independent of (uncertain) distances. 

For R Leo, essentially a single spatial frequency was sampled at H and K 7 bands, and 
the visibilities are sufficiently low that mathematically, three solutions to the UD fit exist: 
one in the main lobe of the visibility function, and two in the second lobe. We therefore use 
apriori knowledge of R Leo's nominal near-infrared size to further constrain our solutions. 
For K'-band, we use the diameter measured by Monnier et al. (2004) (30.3 ± 0.2 mas) as a 
starting value for the fit. For H-band, we use the unpublished value obtained by one of us 
(J. D. Monnier) from aperture masking at the Keck telescope (~ 30 mas), to obtain the two 
first lobe solutions indicated in Table 5. Without spatial frequency coverage that follows the 
shape of the visibility function, it is not possible to discriminate between these solutions, 
and therefore we have droped this star from the analysis that follows. 

As is apparent in Table 5, the basic result of this paper is that the near-infrared diameters 
are found to follow the trend: J diameter < H diameter < K' diameter. This result is also 
illustrated in Figure 2, where the angular diameter ratios are plotted as a function of pulsation 
phase. The mean values of the diameter ratios are: 



2 For each observation we have estimated the pulsation phase at the observation epoch as follows. The 
AFOEV (http://cdsweb.u-strasbg.fr/afoev/) and AAVSO (http://www.aavso.org/) databases were 
used to obtain the epochs of maximum light before and after the IOTA observation, o max i and o maX 2 
respectively. The phase at the IOTA epoch is then estimated as Phase = (IOTA epoch — o max i) / (o maX 2 — 

°ma X 1 ) • 
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R J/H = 0.93 ±0.02 
= 0.89 ±0.02 

(1) 

where the I — a uncertainty is given by the error in the mean, and the (sample RMS, number 
of observations) are (0.09,16) and (0.09,24) for J/H and H/K' respectively. 

We finally comment on a few cases for which there are intriguing changes in the measured 
sizes, above the levels of calibration accuracy demonstrated in § 3. As can be seen in Table 5, 
whenever there is both a significant time difference between two epochs of observation and 
the baselines used were different, apparent size changes in the range 10 — 30% can occur 
(R Aur, U Ori, R Lmi, S Crb, R Ser). Given the variety of possible effects at play however 
(physical changes in the stars, inadequacy of the UD model, stellar asymmetries probed 
for different baseline orientations), it is difficult to assess whether these effects are real or 
arise from unusually high calibration errors. These issues might be resolved with further 
observations. 



5. Discussion 

We have measured the diameters of 23 Mira stars simultaneously in the near-continuum 
near-infrared J, H and K' bands, a first for each star in the sample. The simultaneity of 
the measurements guarantees that the same physical state of the star is measured in each 
band. We find the apparent size relation: J diameter < H diameter < K' diameter, in 
agreement with expectations that the shorter near-infrared filter (J-band) better samples 
the true photospheric continuum, and that the H and K' filters contain progressively more 
opacity contamination from molecules in higher layers (see e.g. Jacob & Scholz 2002). In 
terms of the magnitude of the size difference, or diameter ratios, we find that the stars appear 
to be as much as 24% smaller at J than at H (mean 8%) and as much as 27% smaller at H 
than at K 7 (mean 11%). These magnitudes are in general agreement with the models of Jacob 
& Scholz (2002) and Ireland, Scholz & Wood (2004); although clearly detailed comparisons 
must be made for specific stars and specific models. 

Changes in the near-infrared sizes of individual Mira stars as a function of pulsation 
phase have been detected interferometrically (Perrin et al. 1999; Thompson 2002; Thompson, 
Creech-Eakman & van Belle 2002). No single star in our sample has sufficient phase coverage 
for such a study. However, we have attempted to detect the pulsation by treating the sample 
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as a single "synthetic" Mira. A priori, the most promising approach appeared to be to explore 
this effect using the size ratios, due to their calibration advantages. However, as is apparent 
in Figure 2, no sinusoidal signature is present. We have also explored the possibility that 
the pulsation was not seen in the above approach due to the fact that instantaneous spectral 
type, rather than pulsation phase, might be the relevant parameter. Thus, we have used 
"phased spectral types" (column 4 of Table 5) from the literature (Terrill 1969; Lockwood & 
Wing 1971; Lockwood 1972) to search for a monotonic relation (also for the "synthetic Mira") 
between the measured apparent sizes and the spectral type measured for the corresponding 
pulsation phase, also without success. 

Forsaking the advantages of distance independence then, we show in Figure 3 our mea- 
sured J/H/K' linear sizes as a function of pulsation phase. As dicussed by van Belle, Thomp- 
son & Creech-Eakman (2002), Hipparcos distances to Mira stars are rather uncertain, and 
therefore we have used distances from the literature based on period-luminosity (PL) rela- 
tions, as indicated in Table 1. We note that we do not attempt to apply the PL relations 
using estimates of the absolute K-band magnitudes based on our own fringe data because 
the instrumentation and observing strategy were never intended for accurate photometry. 
Although the data errors are relatively large, and systematic deviations are apparent, the 
data suggests a sinusoidal signature with maximum sizes near phase = 0.5 (minimum light), 
as expected (as the photospheric size increases and cools, the star becomes fainter). The 
statistical significance of the sinusoidal signature is however marginal, with reduced x 2 °f 
0.4, 0.8 and 0.9 at J/H/K' respectively, compared with 1.9, 2.0 and 2.2 for models consisting 
of the constant average diameters. We note that while the observed sinusoid amplitudes 
(50 — 60% peak-to-peak) are comparable to those obtained from model predictions (Ireland, 
Scholz & Wood 2004), they are considerably larger than observed (9 — 26% peak-to-peak) 
when individual Miras are followed through their pulsation cycles (Thompson 2002; Thomp- 
son, Creech-Eakman & van Belle 2002). Also, it can be seen in Figure 3, that the linear 
sizes vary rather in unison, explaining why the pulsation was not apparent in the diameter 
ratios of Figure 2. 

Ultimately, we believe that the failure to more convincingly detect the pulsation using 
the approaches just described is due to the fundamental inadequacy of those methods, as 
several physical mechanisms may contribute dispersion at levels comparable to the expected 
signatures: (a) intrinsic star-to-star differences, (b) our filters are broadband, and therefore 
probe a complex combination of emission from multiple layers of molecular gas just above 
the photosphere, with opacities which also change during the pulsation, making the true 
expected signature very complicated, and in particular potentially masking the photospheric 
layer motions (Ireland, Scholz & Wood 2004), and (c) some theoretical models predict cycle- 
to-cycle non-repeatibility of the photospheric diameters, at significant levels (tens of %; Jacob 
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& Scholz 2002; Ireland, Scholz & Wood 2004), which if present would introduce that much 
dispersion in our apparent size vs. phase relations. On that last point, we note however 
that for the only published multi-cycle measurements to date (Thompson 2002; Thompson, 
Creech-Eakman & van Belle 2002) this effect was not observed; although it appears to be 
present (at 5 — 15% levels) in yet un-published data by the same group (Thompson R. R. 
2004, private communication). In the literature, the success of the "synthetic" Mira method 
is mixed: while van Belle et al. (1996) saw a 20% peak-to-peak change in the effective 
temperatures of a sample of 18 Miras, the effect was not explored in the 22 star sample of 
their follow-up paper van Belle, Thompson & Creech-Eakman (2002), and appears not to 
be present by straighforward plotting of their published numbers. We finally note that the 
best case in our sample in which to search for a pulsation signature of an individual star 
(U Per, period = 320 days) was observed at two epochs almost exactly one cycle apart, and 
the measured UD sizes differ by 10, 12, and 6% at J, H and K 7 respectively, a marginally 
significant size change. 

Several optical interferometers around the world are now equipped to measure accurate 
visibilities in broad or narrow spectral bands, and together span the visible to mid-infrared 
spectrum. Although recent efforts over more limited regions of the spectrum have clearly 
illustrated the power of multi-wavelength spatially resolved observations (Perrin et al. 2004; 
Jacob et al. 2004; Ohnaka 2004; Weiner 2004), global modelling fitting has however not 
been attempted to date, but represents the necessary next step toward a better theoretical 
understanding of Mira star atmospheres and dynamics. 
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Table 1. Basic properties of Mira star sample. 
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Table 1 — Continued 



Name 


RA 


DEC V m j n "V r max 


J 


H K 


Spectral 


Distance 


Distance 




(J2000) 


(J2000) 






Type 


(PC) 


reference 


R Aqr 


23 43 49.5 


-15 17 04.2 6.5 - 11 


-0.1 


-1.1 -1.6 


M7 


272 ± 68 


1 



Note. — V magnitudes and spectral types are from the AAVSO and SIMBAD databases, near-infrared 
magnitudes are from the 2MASS database; as these stars are by definition variable, this information is 
meant to be merely representative. 

References. - (1) van Belle, Thompson & Creech-Eakman (2002); (2) van Belle et al. (1997); (3) 
Young (1995); (4) Jura & Kleinmann (1992); (5) Wyatt & Cahn (1983) 
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Table 2. Log of observations. 
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S Her 
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J/H 


HR6542, HR6065 


Ry Lyr 


98/06/15 


0980.0 


S15N35 


H/K' 


HR7192 


R Aql 


98/06/13 


0978.0 


S15N35 


H/K' 


HR7557 
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Table 2 — Continued 



Name 


Date 


Julian Date 


Baseline 


Filters 


Calibrator(s) 




(yy/mm/dd) 


(-2450000) 








R Peg 


98/10/01 


1088.0 


S15N35 


J/H/K' 


HR8608 


W Peg 


98/10/01 


1088.0 


S15N35 


H/K' 


HR9051 


S Peg 


98/09/30 


1087.0 


S15N35 


J/H/K' 


HR8916 


R Aqr 


97/10/17 


0739.0 


S15N15 


J/H/K' 


HR8980 



Note. — The baseline designations are as follows: SxNy refers to South 
telescope being at the x meter station and the North telescope at the 
y meter station. The [length, azimuth] of the S15N15 and S15N35 baselines 
are [21 m,— 3.8 deg] and [38 m,+18 deg] respectively. 
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Table 3. Calibrator estimated diameters and uncertainties. 



Calibrator Calibrator Spectral Type Angular Diameter Reference 



(HR number) 


(HD number) 




(mas) 




HR82 


HD1671 


F5III 


0.65 ± 0.08 


la 


HR175 


HD3817 


G8III 


1.19 ± 0.15 


3 


HR353 


HD7147 


K4III 


1.59 ± 0.20 


3 


HR470 


HD10110 


K5III 


1.84 ± 0.02 


2 


HR787 


HD16735 


KOII 


1.29 ± 0.16 


3 


HR876 


HD18339 


K3III 


1.45 ± 0.18 


3 


HR978 


HD20277 


G8IV 


1.14 ± 0.09 


4 


HR1133 


HD23193 


A2m 


0.23 ± 0.04 


4 


HR1588 


HD31579 


K4III 


1.78 ± 0.23 


3 


HR1698 


HD33856 


K3III 


2.13 ± 0.02 


lb 


HR1830 


HD36134 


K1III 


1.29 ± 0.16 


3 


HR1866 


HD36678 


MOIII 


2.55 ± 0.03 


2 


HR2047 


HD39587 


GOV 


1.11 ± 0.20 


4 


HR2105 


HD40486 


KO 


1.87 ± 1.33 


4 


HR2169 


HD42049 


K4III 


2.17 ± 0.28 


3 


HR2188 


HD42466 


K1III 


0.96 ± 0.12 


3 


HR2338 


HD45466 


K4III 


2.17 ± 0.28 


3 


HR2864 


HD59294 


K1III 


2.24 ± 0.03 


2 


HR3115 


HD65522 


K2 


1.50 ± 0.98 


4 


HR3319 


HD71250 


M3III 


4.61 ± 0.60 


3 


HR3791 


HD82522 


K4III 


1.13 ± 0.14 


3 


HR3877 


HD84561 


K4III 


2.15 ± 0.27 


3 


HR4081 


HD90040 


K1III 


1.48 ± 0.19 


3 


HR4958 


HD114149 


KOIII 


1.51 ± 0.19 


3 


HR5020 


HD115659 


G8III 


3.27 ± 0.41 


3 


HR5228 


HD121156 


K2III 


0.96 ± 0.12 


3 


HR5312 


HD124206 


K3III 


1.78 ± 0.22 


3 


HR5674 


HD 135438 


K5 


2.57 ± 2.70 


4 


HR5824 


HD139663 


K3III 


1.99 ± 0.02 


2 


HR5940 


HD142980 


K1IV 


1.27 ± 0.16 


3 
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Table 3 — Continued 



Calibrator 


Calibrator 


Spectral Type 


Angular Diameter 


Reference 


(HR number) 


(HD number) 




(mas) 




HR6065 


HD146388 


K3III 


1.22 


± 


0.15 


3 


HR6542 


HD159353 


KOIII 


1.10 


± 


0.14 


3 


HR6770 


HD165760 


G8III 


1.74 


± 


0.22 


3 


HR6793 


HD166229 


K2III 


1.56 


± 


0.20 


3 


HR7192 


HD176670 


K3III 


2.33 


± 


0.03 


2 


HR7557 


HD187642 


A7V 


3.25 


± 


0.41 


lc 


HR7804 


HD194258 


M5III 


6.62 


± 


0.87 


3 


HR8608 


HD214298 


K5 


1.77 


± 


0.87 


4 


HR8916 


HD220954 


K1III 


1.94 


± 


0.02 


2 


HR8961 


HD222107 


G8III 


2.64 


± 


0.33 


3 


HR8980 


HD222547 


K4III 


2.73 


± 


0.35 


3 


HR8987 


HD222643 


K3III 


2.00 


± 


0.25 


3 


HR9010 


HD223173 


K3II 


2.48 


± 


0.32 


3 


HR9029 


HD223559 


K4III 


2.15 


± 


0.27 


3 


HR9051 


HD224128 


K5 


1.49 


± 


0.55 


4 



References. - (1) from CHARM catalog (Richichi & Percheron 2002); (la) 
Lane, Boden & Kulkarni (2001); (lb) Cohen et al. (1999); (lc) van Belle et al. 
(2001); (2) Borde et al. (2002); (3) van Belle (1999); (4) getCal/Fbol 
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Table 4. Calibrated visibility data. 



Name Date Julian Date Filter u v Visibility 

(yymmdd) (arcsec^ 1 ) (arcsec^ 1 ) Amplitude 



R And 


971121 


2450774 


H 


-2.22 


112.25 


0. 


.463 


± 





.035 


R And 


971121 


2450774 


H 


-4.79 


112.16 


0. 


.419 


± 





.035 


R And 


971121 


2450774 


H 


-12.91 


111.48 


0. 


.485 


± 





.035 


R And 


971121 


2450774 


K' 


-13.46 


83.05 


0. 


.766 


± 





.040 


R And 


971121 


2450774 


K' 


-15.92 


82.57 


0. 


.712 


± 


0. 


.040 


Z Cet 


981104 


2451122 


J 


57.79 


118.31 


0. 


.858 


± 





,026 


Z Cet 


981104 


2451122 


J 


54.73 


118.25 


0. 


.850 


± 





.026 


Z Cet 


981104 


2451122 


J 


51.98 


118.19 


0. 


.863 


± 





,026 


Z Cet 


981104 


2451122 


H 


52.04 


89.86 


0. 


.858 


± 





.021 


Z Cet 


981104 


2451122 


H 


50.23 


89.80 


0. 


.843 


± 





,021 


Z Cet 


981104 


2451122 


H 


48.24 


89.74 


0. 


.838 


± 





.021 


Z Cet 


981104 


2451122 


K' 


26.08 


67.09 


0. 


.913 


± 


0. 


.010 


Z Cet 


981104 


2451122 


K' 


24.98 


67.07 


0. 


.918 


± 





.010 


U Per 


971119 


2450772 


J 


-10.34 


140.05 


0. 


.509 


± 





.022 


U Per 


971119 


2450772 


J 


-14.74 


139.53 


0. 


.515 


± 





,022 



Note. — The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. The complete data set is also available in 
the OI-FITS format (Pauls et al. 2004) upon request. 



Table 5. Near- infrared angular diameters and diameter ratios. 



IN dlllt; 


11 i on lof n 

j uiian uaie 




LJ h o con 


D T 


n„ 


UK 


-Kj/H 


Rh/k 








opecciai lype 

1 \/l fill K n n nri 1 

^ivi suociassj 


^mas ) 


(mas) 


(mas) 






R And 


0774.0 


0.6 


6.0 




6.6 ± 0.2 


5.7 ± 0.3 




1.16 ± 0.07 


Z Cet 


1122.0 


0.25 


5.0 


2.7 ± 0.2 


3.5 ± 0.2 


3.7 ± 0.2 


0.76 ± 0.05 


0.95 ± 0.06 


U Per 


0772.0 


0.8 


6.0 


5.0 ± 0.1 


5.4 ± 0.1 


6.6 ± 0.2 


0.93 ± 0.02 


0.83 ± 0.03 


U Per 


1088.0 


0.9 


5.5 


4.5 ± 0.1 


4.85 ± 0.05 


5.0 ± 0.3 


0.93 ± 0.01 


0.97 ± 0.05 


R Per 


1089.0 


0.3 


6.3 


< 1.6 


< 0.9 


< 2.7 






R Per 


1123.0 


0.5 


7.5 


< 1.5 


< 2.2 


< 1.8 






R Aur 


0734.0 


0.7 


8.5 


11.1 ± 0.1 


10.5 ± 0.3 




1.05 ± 0.03 




R Aur 


0772.0 


0.8 


8.5 




9.3 ± 0.1 


10.1 ± 0.1 




0.92 ± 0.01 


R Aur 


0775.0 


0.8 


8.5 




9.2 ± 0.1 


10.1 ± 0.1 




0.91 ± 0.01 


S Ori 


1123.0 


0.1 


6.3 




9.1 ± 0.1 


9.6 ± 0.2 


... 


0.95 ± 0.02 


U Ori 


0739.0 


0.9 


9.0 


10.9 ± 0.1 


12.5 ± 0.1 




0.87 ± 0.01 




U Ori 


0774.0 


0.04 


8.2 


11.5 ± 0.1 










U Ori 


0776.0 


0.04 


8.2 




9.2 ± 0.3 


10.5 ± 0.4 




0.87 ± 0.05 


X Aur 


0736.0 


0.1 


3.4 


< 3.3 


< 3.2 


< 5.2 






R Cnc 


0875.0 


0.4 


7.5 




15.3 ± 0.2 


18.9 ± 0.2 




0.81 ± 0.01 


R LMi 


0773.0 


0.2 


8.0 






12.1 ± 0.1 






R LMi 


0873.0 


0.5 


9.0 


12.5 ± 0.1 


12.8 ± 0.2 


14.2 ± 0.2 


0.97 ± 0.01 


0.91 ± 0.02 


R LMi 


0875.0 


0.5 


9.0 


12.6 ± 0.1 


13.5 ± 0.1 


14.7 ± 0.1 


0.94 ± 0.01 


0.91 ± 0.01 


R Leo 


0874.0 


0.4 


7.5 




23.8 ± 0.3 
32.4 ± 0.4 


29.91 ± 0.27 




0.79 ± 0.01 


R Hya 


0873.0 


0.8 


8.0 


19.8 ± 1.1 


22.7 ± 0.6 


23.9 ± 0.5 


0.87 ± 0.05 


0.95 ± 0.03 



Table 5 — Continued 



Name 


Julian Date 


Phase 


Phased 




D H 

n 


D K 


Rj/H 


It IT 1 1< 




(-245000) 




Spectral Type 
(M subclass) 


(mas) 


(mas) 


(mas) 








R Hya 


0875.0 


0.8 


8.0 


20.5 ± 0.2 


22.0 ± 0.4 


25.8 ± 0.2 


0.93 


± 0.02 


0.86 ± 0.01 


W Hya 


0874.0 


0.6 


7.0 




31.2 ± 0.3 


39.9 ± 0.2 






0.78 ± 0.01 


S Crb 


0874.0 


0.4 


7.5 


9.2 ± 0.7 


8.9 ± 0.9 


9.2 ± 0.9 


1.04 


± 0.13 


0.96 ± 0.14 


S Crb 


0977.0 


0.6 


9.0 




9.1 ± 0.5 


11.2 ± 0.4 






0.81 ± 0.05 


RS Lib 


0873.0 


1.0 


7.5 


9.8 ± 0.4 


9.1 ± 0.9 


9.3 ± 1.1 


1.08 


± 0.11 


0.98 ± 0.15 


R Ser 


0875.0 


0.0 


4.0 


7.5 ± 0.1 


8.3 ± 0.7 


10.4 ± 0.4 


0.90 


± 0.07 


0.80 ± 0.07 


R Ser 


0880.0 


0.1 


6.0 




6.8 ± 0.1 


7.9 ± 0.1 






0.86 ± 0.01 


S Her 


0980.0 


0.6 


8.0 


5.4 ± 0.3 


5.6 ± 0.2 




0.97 


± 0.06 




Ry Lyr 


0980.0 


0.35 


7.7 




2.7 ± 0.3 


< 2.03 








R Aql 


0978.0 


0.9 


7.5 




9.3 ± 0.1 


10.6 ± 0.1 






0.88 ± 0.01 


R Peg 


1088.0 


0.0 


7.0 


6.0 ± 0.4 


7.0 ± 0.3 


9.6 ± 0.2 


0.86 


± 0.06 


0.73 ± 0.03 


W Peg 


1088.0 


0.1 


6.9 




8.2 ± 0.1 


9.6 ± 0.1 






0.86 ± 0.01 


SPeg 


1087.0 


0.1 


5.2 


3.9 ± 0.1 


4.9 ± 0.3 


4.8 ± 0.5 


0.79 


± 0.05 


1.02 ± 0.12 


R Aqr 


0739.0 


0.4 


7.5 


17.7 ± 0.2 


17.7 ± 0.1 


20.8 ± 0.8 


1.00 


± 0.01 


0.85 ± 0.03 



Note. - For R Leo, the two H-band diameter values indicated correspond to the two possible solutions that 
exist in the second lobe of the visibility function, see § 4 for discussion. 
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Fig. 1. — Example visibility data. The four panels illustrate the general finding that lower 
visibilities are measured at the longer near-infrared wavelengths (J diameter < H diameter 
< K' diameter). The bottom-right panel also illustrates the data precision. In all panels, the 
abscissa is the spatial frequency given by the projected baseline B v j Ao = \Ju 2 + v 2 ; and the 
ordinate is the visibility amplitude (1.0 for a point-like source, 0.0 for a completely resolved 
extended source). 
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Fig. 2. — Angular diameter ratios. For each panel, the solid line indicates the measured 
mean ratio, the dashed lines indicate the error in the mean, and the dotted lines indicate 
the RMS of the sample. 
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Fig. 3. — Linear diameters as a function of pulsation phase. In each panel, the dashed 
line indicates a sinusoidal fit to the data. The sinusoidal fits have 4 free parameters: mean, 
amplitude, period and phase; except for the K'-band data, where we fix the period and phase 
to the values found in the J and H fits. See 5 5 for discussion. 



